y>A /~37 
/t 10 2 - 

NASA Army Research Laboratory 

Technical Memorandum 106678 Memorandum Report ARL-TR-520 

A /x 


Analysis of the Effects of Surface Pitting 


and Wear on the Vibrations of a 
Gear Transmission System 

F.K. Choy and V. Polyshchuk 
The University of Akron 
Akron, Ohio 

J.J. Zakrajsek 
Lewis Research Center 
Cleveland, Ohio 

R.F. Handschuh 
Vehicle Propulsion Directorate 
U.S. Army Research Laboratory 
Lewis Research Center 
Cleveland, Ohio 

D.P. Townsend 
Lewis Research Center 
Cleveland, Ohio 


Prepared for 
AUSTRIB ’94 

sponsored by the Institution of Engineers, Australia 
Perth, Western Australia, December 5-8, 1994 


rH 


CsJ 


s O 


<\| 

to 


fTJ 

1 

f— 

un 

U 

o 

c 

z 

Z) 



a: 




< 



UJ 

UJ 



X 

3: 



b- 


1 



O 

<v 


u. 

Z 

-j 


o 

< 




QC 

• 


to 

O < 

< 



2! UU 

CO 


CO 

o 

< 


> 

b- 



-J 

b- < 

w 

Cl 

< 

M 




a cl 


O 

< 

a 

X 

H 


UJ 

UJ 



u to 

b- 



< 2! 

to 


CO 

CL O 

> 

U 

r- 

a: ^ 

in 

QJ 

>0 

ID ►— 


4J 

o 

CO < 

Z 

C 

o 

QC 

o 

0) 

rH 

IL CD 

►— « 

u 

1 

a ~ 

00 


s: 

> 

00 

x: 

h- 

00 

»— 4 

o 


I b- liJ I L 

< U I 13 

W UU h Z 0) 

< ul < i n 

z Ik z o: CL 

^ oj O h it 


Csj 

vO 

r- 

o 

t-4 

o 

o 


N 

m 


m 

O 



National Aeronautics and 

Space Administration RESEARCH LABORATORY 



U.S. ARMY 


Analysis of the Effects of Surface Pitting and Wear on the 
Vibrations of a Gear Transmission System 

F. K. Choy*,V. Polyshchuk*, J J. Zakrajsek^, R.F. Handschuh?, and D.P. Townsend^ 

♦The University of Akron, tNational Aeronautics and Space Administration, 
t Vehicle Propulsion Directorate, U.S. Army Research Laboratory, Lewis Research Center 


Abstract 

A comprehensive procedure to simulate and analyze the vibrations in a gear transmission system with surface pitting, wear, 
and partial tooth fracture of the gear teeth is presented. An analytical model was developed where the effects of surface 
pitting and wear of the gear tooth were simulated by phase and magnitude changes in the gear mesh stiffness. Changes in the 
gear mesh stiffness were incorporated into each gear-shaft model during the global dynamic simulation of the system. The 
overall dynamics of the system were evaluated by solving for the transient dynamics of each shaft system simultaneously 
with the vibration of the gearbox structure. In order to reduce the number of degrees-of-ffeedom in the system, a modal 
synthesis procedure was used in the global transient dynamic analysis of the overall transmission system. An FFT procedure 
was used to transform the averaged time signal into the frequency domain for signature analysis. In addition, the Wigner- 
Ville distribution was also introduced to examine the gear vibration in the joint time-frequency domain for vibration pattern 
recognition. Experimental results obtained from a gear fatigue test rig at NASA Lewis Research Center were used to evalu- 


1. Introduction 

In the last two decades, problems arising from excessive 
wear and gear tooth surface pitting in gear transmission 
systems have been of increasing concern for a variety of 
gear users. At increased power and higher speeds, gear wear 
and fatigue failures are major concerns in aerospace appli- 
cations. Although regular visual inspections and preven- 
tive maintenance can help to reduce the failure rate of gear 
systems, the cost and down time required make such pro- 
grams inefficient and uneconomical. 

Vibration signature analysis methodologies are being 
developed to non-intrusively examine the health and wear 
of gear transmission systems. Using spectral analysis, the 
amplitude of the frequency spectrum of the measured 
vibration signal is calculated and displayed in a continuous 
manner. However, the spectral analysis technique is diffi- 
cult to apply in a highly complex system where the large 
number of spectral lines often makes it difficult to detect 
significant changes in the spectrum. Another methodology 
is the joint time-frequency approach which applies the 
Wigner-Ville distribution (WVD) [1-3] on the time vibra- 
tion signal of the system. Unlike the regular Fourier trans- 
form process, the WVD provides an instantaneous fre- 
quency spectrum of the system at any instant throughout 
the sampling period (while FFT provides a averaged fre- 
quency spectrum of the total sampling period). The spec- 
tral density of the fundamental exciting frequency and its 
sidebands change as the shaft rotates through a complete 
revolution. Some success has been achieved in applying 
the WVD concept in the health monitoring of gear trans- 
mission systems [3-6]. However, a complete vibration sig- 
nature database is needed for development of an effective 
pattern recognition scheme. In order to populate such data- 
bases, the development of an accurate analytical procedure 
to predict vibrations in gear systems due to wear and fatigue 
failure is necessary. 


The objective of this paper is to develop a comprehensive 
procedure to simulate and analyze the vibration in a gear 
transmission system with effects of surface pitting and wear 
of the gear teeth under normal operating conditions. The 
effects of changes in magnitude and phase of the mesh stiff- 
ness at one particular tooth or a number of consecutive teeth 
were evaluated in order to simulate the effects of surface 
pitting and wear. The effects of these localized changes in 
the gear mesh were incorporated into each gear-rotor model 
for the dynamic simulation [7-9]. The dynamics of each 
gear-rotor system were coupled with each other through 
the gear mesh interacting forces. The coupling between the 
rotors and the casing structure were generated through the 
bearing support forces. The global vibrations of the system 
were evaluated by solving the transient dynamics of each 
rotor system simultaneously with the vibration of the cas- 
ing. In order to minimize the computational effort, the num- 
ber of degrees-of-freedom of the system were reduced by 
using a modal synthesis procedure [7,8]. Experimental 
vibration results obtained from a gear fatigue test rig at 
NASA Lewis Research Center [11] were used to verify the 
analytical procedure. 

2. Solution Procedures 

2.1 Dynamics of the Gear-Shaft Configuration and the 
Gearbox System 

The dynamics of the ith individual gear-shaft system can 
be evaluated through the equations of motion for the vibra- 
tions of a individual rotor-bearing-gear system as shown in 
Figure 1 [7,8], given in matrix form, as 

[M]{w j } + [K s J{w i } = {F bi (t)} + {F gi (t)} + {F uj (t)} (1) 

where [M] and [Kg] are respectively the mass and shaft stiff- 
ness matrices of the rotor, {W-} is the general displacement 
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vector of the ith rotor in the its local coordinate system, 
and, {F bi (t)}, {F gi (t)}, and {F ui (t)} are respectively the 
force vectors acting on the ith rotor system due to bearing 
forces, gear mesh interactions, and mass-imbalances. 

The equations of motion of the gearbox with p rotor sys- 
tems can be expressed as 

K] {wj + M { w c}= ( 2 ) 

i=l 

where [T c * ] represents the coordinate transformation 
between the ith rotor and the gearbox. 

2.2 Evaluation of Bearing Forces 

The bearing forces {F^-(t)} for the ith rotor can be evalu- 
ated as 

taWj-Mtw.HTjWa}) 

WW-WW ® 

where [C^-] and [K^j] are respectively the damping and 
stiffness of the bearing, [T- c ] is the coordinate transforma- 
tion matrix for the gearbox with respect to the ith rotor, and 
W c - are the casing displacements at the rotor locations. 

2.3 Evaluation of Gear Forces 

The gear forces generated from the gear mesh interac- 
tion^] can be written as 

< 4 > 

where {F^(t)} is the vector containing the gear forces and 
moments resulting from the relative rotation between the 
two mating gears and {F^(t)} is the vector containing gear 
forces and moments due to the translational motion between 
the two gears. 

2.4 Modal Synthesis Procedure 

In order to calculate the transient and steady state dynam- 
ics of the system, the coupled rotor and casing equations of 
motion must be solved simultaneously. To minimize the com- 
putational effort, the modal transformation[7,8] procedure 
will be applied to reduce the degrees of freedom of the 
global equations of motion. Using m undamped mode 
shapes of the ith rotor system [fyj, <J)|2* ^3* — * ^ 

m c undamped mode shapes of the gearbox $ 03 ’ 

<t> cmc ], the rotor displacement for the ith rotor can be 
written as 

m 

m-IaM <» 

j=l 


and, similarly, the gearbox displacements as 

{W c H*c]{A c } ( 6 ) 

where { Aj} and { A c } are the modal time functions of the 
ith rotor and the gearbox respectively. Using the expansion 
in equation (5), the equations of motion for the ith rotor in 
equation (1) can be written as 

[m] [<t> c ] {Aj}+[K S ] [<^] {Aj} 

={ F bi«Myt)} + {Fjt)} (?) 

Premultiplying by [<(>j] T and using the orthogonality condi- 
tions of the mode shapes [7], the ith rotor equations of 
motion can be written as 

fthMfAiHsy+foM?*} <8> 

where [ A^] is the diagonal matrix of the squares of the 
natural frequencies of the system. 

For the gearbox system, a similar transformation is carried 
out as equation (2) can be written as 

(aJ+KIKH'W » 

For a system of k rotors, equation (8) can be repeated k 
times and solved with the casing equation (9) simulta- 
neously for the modal accelerations {Aj} and { A c }. A 
numerical integration scheme is used to integrate the 
accelerations to obtain velocities and displacements at each 
time step for transient calculations [7]. 

3. Signature Analysis of Vibration Signal 

3.1 Frequency Domain Analysis 

The frequency spectrum analysis is used by applying a dis- 
crete Fourier Transform on the average time signal x(t) such 
that the spectral components are 

X(k) = T£x(t)exp[^^] (10) 

i=0 

where x(t) is the time averaged of the vibration signal W(t) 
and T is the sampling interval. The frequency components 
are examined in the frequency domain and compared with 
those obtained at various stages of the fault development 
in the experimental gear test rig. 

3.2 Joint Time-Frequency Analysis : The Wigner-Ville 
Distribution 

To examine the vibration signal in a joint time-frequency 
domain, the Wigner-Ville method [1-3] is used in this study. 
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The Wigner-Ville distribution will provide an inter-domain 
relationship between time and frequency during the period 
of the time data window. The WVD (Wigner-Ville Distri- 
bution) can be written as: 

oo 

WV(t, f ) = J x^t + 1 j x * ^t - 1 J e" j2;rft dt (1 1) 

— oo 

where WV(t, f) is the Wigner-Ville distribution in both the 
time domain t and the frequency domain f. To allow sam- 
pling at the Nyquist rate and eliminate the concentration of 
energy around the frequency origin due to the cross prod- 
uct between negative and positive frequency [1,2], the ana- 
lytic signal was used in evaluating the WVD. The analytic 
signal s(t) is defined as 

S(t) = x(t) + jH[x(t)] (12) 

Where H[x(t)] is the Hilbert transform of x(t). 

In order to avoid a repetition in the time domain WVD, a 
weighting function[4] is added to the time data before the 
evaluation process. Such process may decrease the resolu- 
tion of the distribution, but it will eliminate the repetition 
of peaks in the time domain, and, thus the interpretation of 
the result will be substantially easier. 

4. Description of Experimental Study 

The experiment was performed on the spiral bevel gear 
fatigue test rig [1 1], as illustrated in Figure 2, at the NASA 
Lewis Research Center. The primary purpose of this rig is 
to study the effects of gear tooth design, gear materials, 
and lubricants on the fatigue strength of aircraft quality 
gears. Because spiral bevel gears are used extensively in 
helicopter transmissions to transfer power between non- 
parallel intersecting shafts, the use of this fatigue rig for 
diagnostic studies is practical. Vibration data from an 
accelerometer mounted on the pinion shaft bearing hous- 
ing was captured by an analog to digital conversion board. 
The 12- tooth test pinion, and the 36-tooth gear have a 35 
degree spiral angle, a 1 in. face with, a 90 degree shaft 
angle, and 22.5 degree pressure angle. The pinion trans- 
mits 720 hp at a nominal speed of 14,400 rpm. The test rig 
was stopped several times during the test for gear damage 
inspection. The test was concluded at 17.8 operational hours 
when a broken tooth was detected visually during one of 
the shutdowns. 

Pictures of tooth damage on the pinion at various stages in 
the test are shown in Figure 3. At the first rig shut-down, at 
about 5.5 hours into the test, a small pit was observed on 
one of the teeth on the test pinion, as illustrated in Figure 
3 A. The test was stopped again at approximately 12 hours 
and the pitted area spread to cover approximately 75% of 
the face of the pinion tooth, as seen in Figure 3B. In addi- 
tion, pitting started to appear on the adjacent teeth. Figure 
3B shows the pinion at the end of the test, at 17.8 hours. It 
was found that one of the three heavily pitted pinion teeth 
had experienced a tooth breakage, losing one third of the 
tooth, as shown in the figure. 


5. Discussions of Results 

To study the effects of gear tooth pitting and wear on the 
dynamics of the rotor system, the numerical simulation 
procedure described above was used to model the vibra- 
tions of the pinion gear in the test rig. During the experi- 
mental study, vertical direction vibration signals from the 
pinion gear are time synchronously averaged for spectral 
analysis and analysis using the joint time-frequency 
distribution(WVD). In order to perform an accurate compari- 
son, the averaged time signal from the vertical vibration of 
the pinion gear is also generated using the numerical model. 
During these simulations, approximate gear mesh stiffness 
models are developed to simulated the effects of wear and 
pitting of the pinion tooth on the dynamics of the system. 

As it has been established, the changes due to gear tooth 
wear or failure can be represented by the amplitude and 
phase changes in vibration, which, in turn, can be repre- 
sented by magnitude and phase changes in mesh stiff- 
ness^]. To demonstrate the effects of mesh stiffness 
change on gear vibration, the variation of the mesh stiff- 
ness model used for this study is given in Fig 4. The 
“undamaged” configuration of the mesh stiffness is given 
by 0 degree phase change(Fig. 4A), and 0 % amplitude 
reduction(Fig. 4B).During the wear and pitting process, two 
types of stiffness changes are examined, i.e., the phase 
changes, shown in Figure 4A, and the amplitude changes, 
shown in Figure 4B. Figures 5 and 6 show the time, fre- 
quency, and joint time-frequency analysis(WVD) of the 
pinion gear vibration signals with the approximated changes 
in gear mesh stiffness. 

Figure 5 shows the effects of mesh stiffness phase changes 
in the WVD representation of the predicted vibration sig- 
nal. As seen in Figure 5, a phase change in the mesh stiff- 
ness at the 6th tooth of the 12-teeth pinion resulted in a 
temporary increase of amplitude and phase of the pinion 
vibration time signal during the 6th tooth pass location. As 
the phase shift in the mesh stiffness increases, from 
1.5 degrees to 4.5 degrees, the changes in amplitude and 
phase in the vibration signal become more pronounced. In 
the frequency spectra, this change in mesh stiffness will 
result in the increase of the amplitude in the sideband 
frequencies. However, as discussed earlier, although the 
frequency spectrum provides good indications of the exist- 
ence of the non-synchronous components, it can not dis- 
tinguish the time locations of their occurrences. The joint 
time-frequency analysis using WVD, shows the existence 
of various frequency components as the pinion rotates 
through a complete revolution of 360 degrees. Note that, 
in this case, the WVD shows a continuous excitation of the 
mesh frequency ( 12 x rotational speed) throughout the com- 
plete 360-degree revolution while subsynchronous com- 
ponents of 8 times, 4 times, and 1 times rotational speed 
are occurring at the 6th and 7th tooth pass locations. 

Figure 6 shows the effects of reductions in mesh stiffness 
at the 6th gear location. With the reduction of mesh stiffness, 
a substantial change in the vibration at the 6th tooth pass 
location (150 - 180 degrees) is observed. Note that in the 
case of 50% stiffness reduction, the time vibration amplitude 
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at the 6th tooth pass location(at approximately the 180 de- 
gree location for the 12 teeth pinion) almost vanishes and a 
much larger amplitude at the 7th tooth pass location is gen- 
erated. In addition, the vibration amplitudes of the 8th and 
9th tooth pass locations are reduced. These reductions in 
vibration amplitudes at mesh frequency resulted in a much 
higher sub-synchronous components in the frequency spec- 
trum as shown in Figure 6C. The WVD shows a distinct 
type of cross pattern at the intersection of the mesh fre- 
quency and the 6th tooth pass location with a continuous 
mesh frequency component throughout the complete pin- 
ion revolution. 

Figure 7 shows the pinion gear vibration signature analysis 
of the experimental time signal acquired at A) 12 hours 
when one tooth is severely pitted (Figure 3B), and B) 17.8 
hours, when three consecutive teeth are pitted and one has 
a tooth fracture (Figure 3C). To numerically simulate these 
phenomena, two gear mesh stiffness models, as shown in 
Figure 8 A and 8B, which include a combination of phase 
shift and amplitude change, are introduced. Figure 8A rep- 
resents the gear mesh stiffness for a heavily pitted tooth 
during a single tooth pass. The mesh stiffness is simulated 
by a 50% loss in stiffness at approximately the first 20% of 
the tooth contacting period. Figure 8B represent the mesh 
stiffness, for a three tooth pass period, consisting of one 
broken tooth with two heavily pitted teeth at the adjacent 
sides. Note that the stiffness of the middle (broken) tooth is 
simulated by a 50% loss of stiffness, while the other 
adjacent(pitted) teeth are simulated by the stiffness reduc- 
tion similar to that of the single tooth case as shown in Fig. 
8A. Additional frictional effects are also added into the 
model to simulate the roughness of the tooth surface due to 
pitting. The simulated vibration signature of the pinion gear 
is given in Figure 9. Comparing Figures 9A and 7a, for the 
single tooth damage case at 12 hours, one may notice the 
similarities between both the frequency spectra and the 
WVD display. Some of the unevenness in the experimen- 
tal time signal is mainly due to the modulation of frequen- 
cies due to other excitations in the test rig which are not 
numerically modeled. For the tooth break-off case at 17.8 
hours, Figures 7B and 9B, both the numerical and the ex- 
perimental WVD display a large cross pattern at the 6th 
tooth pass location due to tooth break-off. However, some 
discrepancies have been detected between the experimen- 
tal and the numerical time signal at the 4th and 5th tooth 
pass locations. The experimental time signal consists of 
some higher frequency, smaller amplitude vibration modu- 
lation, which are not being numerically simulated. This 
additional modulated signal resulted in the excitation of 
the 14 times rotational speed component, as shown in the 
frequency spectrum in Figure 7B, and, also, in turn, is 
responsible for the small differences created in the WVD. 

6. Summary and Conclusions 

A numerical procedure has been developed to simulate the 
dynamics of gear transmission systems with the effects of 
gear tooth damage due to wear and pitting. The work pre- 
sented in this paper can be summarized as follows: 


1) A modal synthesis methodology has been developed to 
simulate the dynamics of gear transmission systems. While 
the computational efforts has been greatly reduced by modal 
transformation, the numerical results generated maintain 
good accuracy. 

2) Gear tooth damage due to wear and pitting can be simu- 
lated by amplitude and phase changes in the gear mesh stiff- 
ness model. The gear mesh model developed can easily be 
incorporated into the global transmission system for 
dynamic predictions. 

3) Using the time averaging technique, frequency spec- 
trum analysis, and the Wigner-Ville distribution, a signa- 
ture analysis scheme can be developed to examine and char- 
acterize the vibration signal of the gear system. 

4) A parametric study of the effects on the vibration signal 
due to various changes in the gear mesh stiffness model, 
simulating various degrees of pitting and wear damage, 
could provide a comprehensive database for gear fault 
detection and damage estimation research. 
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Figure 1 . — Schematic of the rotor-gear bearing system. 


Figure 2. — Picture of the bevel gear test rig in NASA Lewis 
Research Center. 





Figure 3. — Pictures of the damaged pinion teeth, (a) 5.5 hr. (b) 12 hr. (c) 17.8 hr. 




Figure 4. — Stiffness changes in the gear mesh model. 
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Figure 5. — Simulated pinion vibration signature due to phase 
change in the gear mesh stiffness, (a) 1.5°. (b) 3.0°. (c) 4.5°. 
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Figure 6. — Simulated pinion vibration signature due to amplitude 
change in the gear mesh stiffness, (a) 0% reduction, (b) 25% 
reduction, (c) 50% reduction. 
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Figure 7. — Experimental pinion vibration signature due to damage on pinion teeth due to wear and pitting, (a) Single tooth, (1 2 hr), 
(b) Three teeth, (17.8 hr). 



Figure 8. — Gear mesh stiffness model 
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Figure 9. — Numerically simulated pinion vibration signature due 
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damage on pinion teeth due to wear and pitting, (a) Single tooth. 
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